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Abstract 

In this paper we consider a hybrid elastic model consisting of a Timoshenko beam and 
a tip load at the free end of the beam. Under the eqnal speed wave propagation condition, 
we show polynomial decay for the model which inclndes the rotary inertia of the tip load 
when feedback bonndary moment and force controls are applied at the point of contact 
between the beam and the tip load. 


1 Introduction 

Beam structures have been studied extensively in the last decades: Euler-Bernoulli, Rayleigh 
and Timoshenko beams. The latest model is more accurate since it takes into account not 
only the rotary initial energy but also its deformation due to shear (see Timoshenko’s book for 
physical explanations: [23]). A non-exhaustive list of contributions is: |3j, jlj, [6], [7], [H], [T2] . 

m, na, na, m, m- 

In this paper, we study the stabilization of a Timoshenko beam which has a tip load attached 
to one free end. The beam is clamped at one end while the tip load is hxed to the other end 
a: = 1 in such a manner that the center of mass of the load is coincident with its point of 
attachment to the beam. We assume interaction between the beam and the load. Thus the 
forces and moments within the vibrating beam are transmitted to the tip load which moves in 
accordance with Newton’s law. Dissipation is introduced into the coupled model by applying 
feedback boundary moment and force controls on the displacement and shear velocities. Mul¬ 
tiplying the initial equations by suitable constants and rescaling in time, the coupled motions 
of the beam-load structure are governed by the following problem : 
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(1) 

{utt-{Ux + y)x){.x,t) = 

0, 

for 

(x,t) e (0,1) 

X (0,oo)^ 

(2) 

{ytt - ayXX + b{ux + y)){x, t) = 

0, 

for 

{x,t) e (0,1) 

X (0,cx))^ 

(3) 

u{0,t) = y{0,t) = 

0, 

for 

t e (0, oo). 


with the initial conditions 





(4) 

u{x,0) = uq{x), ut{x,0) = ui{x), y{x,{ 

3) = 

= yo{x), yt{x,^) = yi{x), 

for X G 

and 

the boundary dissipation law 





(5) 

t) + ki{ux{l, t) + y{l, t)) ■■ 


-k2Ut{l,t), 

for t G (0, ( 

X)), 

(6) 

ytt{lE) + hyx{lE) 

= - 

-kiyt{l,t), 

for t G (0, 

oo), 


where a, b, ki, k 2 , ks, k^ are strictly positive constants. 

Denote by p, Jp, El, n, u{x,t) and ip{x,t), the mass density, the moment of mass inertia, the 
rigidity coefficient, the shear modulus of the elastic beam, the lateral displacement at location 
X and time t and the bending angle at location x and time t respectively. Then, our model co¬ 
incides with those of 0,0, HH, EU, ... with u{x, t) = u i X, 


a = and b = 

nE E 


-t ) , y{x,t) = -ip 


X, 


^ t 


This system is studied by Kim and Renardy (|T2]), but with other boundary dissipation laws 
and it is then proved to be exponentially stable. 

M. Bassam, D. Mercier, S. Nicaise and A. Wehbe also consider the same system but with other 
boundary dissipation laws. They study the decay rate of the energy of the Timoshenko beam 
with one boundary control acting in the rot at ion-angle equation. Under the equal speed wave 
propagation condition (a = 1) and if b is outside a discrete set of exceptional values, using 
a spectral analysis, the authors prove non-uniform stability and obtain the optimal polyno¬ 
mial energy decay rate. On the other hand, if ^/a is a rational number and if b is outside 
another discrete set of exceptional values, they also show a polynomial-type decay rate using 
a frequency domain approach. See |1] and the references therein, particularly papers by F. 
Alabau-Boussouira (E)), J.E. Munoz Rivera and R. Racke, papers by S.A. Messaoudi and M.I. 
Mustafa, papers by A. Wehbe and his co-authors: A. Soufyane and W. Youssef... 

The stabilization of the Timoshenko beam is a subject of interest for many other authors re¬ 
cently: D. Feng, W. Zhang with a nonlinear feedback control (0), W. He, S. Zhang, S. Ge (see 
m), O. Morgiil with a dynamic boundary control (|18)). 

The spectral analysis is studied by M.A. Shubov (El] and Q.P. Vu, J.M. Wang, G.Q. Xu, S.P. 
Yung (ED). 

Systems of Timoshenko beams, serially connected or forming a tree-shaped network are another 
interesting point: see Hi. H, BSI, m- 
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The system we consider is also studied by M. Grobbelaar-Van Dalsen in [8] with the same 
feedback controls as ours. It is proved that uniform stability holds under a condition (called 
condition Z.) Unfortunately this condition is not easy to check and the exponential stability 
(for a = 1) remains an open question. This is why, in the present work, we consider the same 
problem which is still open. The main goal of this paper is to prove that the decay of the energy 
is not exponential, but polynomial. 

We conjecture that the same results hold in the case a ^ 1. The computations are more com¬ 
plicated and still have to be performed. 


In Section]^ the abstract framework is introduced and the operator is proved to be m-dissipative 
in the energy space. The existence and uniqueness of a solution of the abstract evolution prob¬ 
lem in appropriate spaces is established. The energy of the solution is then proved to decay to 
zero, using Benchimol Theorem (j3]) (i.e. the operator is proved to have no purely imaginary 
eigenvalues). 

Section is dedicated to a thorough analysis of the spectrum of both the dissipative operator 
and the conservative associated operator. In particular, we give asymptotic expansions for the 
eigenvalues (cf. (36), (37), (38) and (39)). 

It is proved, in Section]^ that the system of generalized eigenvectors of the dissipative operator 
(introduced in the latest section) forms a Riesz basis of the energy space. To this end, we use 
Theorem 1.2.10 of [1] which is a rewriting of Guo’s version of Bari Theorem with another proof 
(see m)- The proof requires the asymptotic analysis performed before. 

At last, the solution is explicitly expressed using the Riesz basis to prove that the energy decays 
polynomially (see Section]^. 

To examplify and validate our results, we give numerical computations and figures representing 
the spectrum of the dissipative operator in Section 


2 Well-posedness and strong stability 

In this section we study the existence, uniqueness and strong stability of the solution of System 
(0-(g. Setting 

G ;= (0,1) and Hl{n) := {/ G H\n) : /(O) = 0}, 
we define the energy space H as follows 

n := Hl{Q) X L\Q) x Hl{Q) x L\Q) x C x C, 


with the inner product defined by 


(7) 


<U,Ui>n := 


vvi + b ^zzi + ab ^i/xVix + + y){uix + yi)^{x)dx 


-j-yyi + ;p77i, 

ki fcs 


for all U = {u,v,y,z,r],'y), Ui = (ui, ui, |/i, zi, 771 , 71 ) G T-i. 

1 1 ^ 

Remark 2.1. The norm <U,U induced by ( ?) is equivalent to the usual norm ofTL. 
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For shortness we denote by || ■ || the L^(r2)-norni. 

Now we dehne the linear unbounded operator A : ^(A) —>■ H by: 

D{A) := {U = {u,v,y,z,r],'y) eU-.u.ye H‘^{n),v G Hl{n),z G HI{VL), 
V = 


yu G D{A), 

(8) AU := (n, {u^ + y)^, z, ay^^ - h{u^ + y), -ki{u^{l) + |/(1)) - k 2 r], -k^y^{l) - k^;^). 

The associated conservative operator is ^ D{A) "H dehned as A but with k 2 = k^ = 0 i.e. 

(9) AU = AqU - k2r]e^ - k^'ye^, 

where U = {u, v, y, z, y, 7 ) G D{A), 65 := (0, 0, 0, 0,1, 0) and eg := (0, 0, 0,0, 0,1). 

System Q-([^ is formally rewritten as the evolution equation 


( 10 ) 


(P) 


Ut{t) = AU{t),t G (0; + 00 ), 
17(0) = f/o, Uo G n, 


with U{t) = {u,ut,y,yt,ut{l),yt{l)) (note that the notation U is kept for this function of the 
time t). 

Proposition 2 . 2 . The operator A is m-dissipative in the energy space TL. 

Proof. We start with the dissipativeness. 

Let U = {u,v,y, z,^,^) G D{A). Using ([^ and (|^, we obtain : 


< AU, U>n= + y)xV + b ^ (^ay^^ - b{u^ + y)^ z 


+ ab ^z^y^ + + z){u^ + y)^ {x)dx 

+ ^(-^i(wa;(l) + 2/(1) - k2y)fj + —(-/C32/x(1) - k4'y)y. 
ki ks 

Then, integrating by parts and using the boundary conditions, we get 


( 11 ) 


h 


ki I 


^<Au,u >H= - ?k(i)r < 0. 


ki 


k?, 


Therefore, A is dissipative. 

Now, we prove that 7l is maximal. For that purpose, we consider any / = (/i, / 2 , /s, fi, 71 , 71 ) 
G PL and we look for a unique element U = {u, v, y, z, rj, 7 ) G D{A) such that 

AU = f. 

Equivalently, we get v = fi, z = fs,rj = /i(l ),7 = / 3 ( 1 ), and we have the following system to 
solve: 

( 12 ) {Ux + y)x{x) = f 2 {x), 

(13) {ayxx - b{ux + 2 /))(a;) = fi{x), 
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(14) 

(15) 


ki{ux{l) + 1 /( 1 )) - k2r] = r]i, 

- k^y^il) - /c 47 = 7i. 


From ( 12 ) it follows Ux{x) +y{x) = F2{x) + ai, where F2{x) = f2{u)du and oi is a constant. 
Consequently ( 13 ) becomes 

( 16 ) ay^^(x) = f4{x) + bF2{x) + hai. 

Let G4 (resp. G2) be the unique solution of {G4)xx = fi (resp. (6*2)11 = -^2) satisfying 
6*4(0) = (6*4 )j:(0) = 0 (resp. 6^2(0) = (G2)a:(0) = 0 ). Then, we hnd that the solutions of ( 16 ) 
satisfying //(O) = 0 are 

. . b x‘^ 

y[x) = -64(0;) + -G2[x) + -oi— + a2X, 
a a a 2 

where 02 is a constant. Now, let y as previously. Clearly y G fl F[\{Vt) and we hnd that 

px 

necessarily u{x) = / {—y{u) + F2{u))du + aix (since m(0) = 0). 


Inserting UxiX) + 2/(1) = -^ 2 ( 1 ) + Oi in (14) we get an equation with only the unknown Oi and 
this equation admits a unique solution. Therefore (15) becomes an equation with a unique 
solution 02 . Finally, inserting these two constants in u and y, it is easy to check that we have 
found a unique U = (u, n, ?/, rj, 7 ) G D{A) such that AU = f. 

Therefore we deduce that 0 G p{A). Then, by the resolvent identity, for A > 0 small enough, 
R{XI — A) = 1-1 (see Theorem 1.2.4 in [H]). □ 


Due to Lumer-Phillips Theorem (see |19( . Theorem 1.4.3), it follows from Proposition 2.2 that 
the operator A generates a C*o-semigroup of contractions on R. Consequently it holds: 

Theorem 2.3. (Existence and uniqueness) 

(1) If Uo & R, then System (P) has a unique solution 

U G G^{R+,R). 

(2) If Uo G D{A), then system (P) has a unique solution 

U G G^{R+,D{A)) nG\R+,R). 

Remark 2.4. Let (Pq) be the conservative problem associated to problem (P) (in other words 


(Po) is Problem (P) with ^2 = ^4 = Oj and Aq he the associated operator then Proposition 2.2 
(resp. Theorem 2.3) remains true for Aq (resp. (Pq) ). 

To end this section we give a hrst stability result: 

Theorem 2.5. (Strong stability) 

System ([^-([^ is strongly stable, i.e for any solution U of (P) with initial data Uq G R, it holds 

lim E{t) =0, 

t^OO 


where E{t) = -||P(t)||^. 
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Proof. Since the resolvent of A is compact in PL, using Benchimol Theorem [5], System (P) is 
strongly stable if and only if A does not have purely imaginary eigenvalues. We have already 
seen that A is invertible. Thus we consider A G M* and U = {u,v,y, z,ri,'y) G D{A) such that 


AU = iXU. 

Since 3? < AU, U >= 0, we get from (0 that T] = n(l) = 0 and 7 = z{l) = 0 , and we deduce 
that (u, v) satisfies 


(17) 

with the boundary conditions 

(18) 


{u^x + yx + X‘^u){x) 


= 0 , 


Now, from the second equation of (17), it follows: 
solution of 


— bux - 

-by + X‘^y){x) = 0, 

u( 0 ) 

= 9(0) 

= 0 , 

n(l) 

= 2 /( 1 ) 

= 0 , 


= w(i) 

= 0 . 

^xx 

{x)—X‘^u 

[x). Thusa|/a;a,(a;) 


aVx 


—au^x\x)—aX'^Ux{x). 
= hux{x) + hy{x) — X^y{x). Then u is 


(19) 


au 


(4)/ 


{x) + (a + l)X^u^^\x) + (A^ — h)X^u{x) = 0. 


Note that, from the boundary conditions (18) and the relations 0 . it also holds u^xi^) = 
UxxxiX = 0- Thus u is solution of (19) and satisfies m(1) = UxiX) = UxxiX) = UxxxiX) = 0. 
Therefore, from the general theory of ordinary differential equations, we deduce that u = 0. 

It follows that y = Q and finally P = 0. Consequently, A has no eigenvalue on the imaginary 
axis. □ 


3 Spectrum analysis for the case a = 1 

3.1 Main results and notation 

Let us begin with announcing the main results concerning the spectrum analysis. The following 
theorem is also a way to introduce the notation which is used during the whole section. That 
is why it is given first whereas establishing its proof is the goal of the following subsections. 

Theorem 3.1. (Spectrum and eigenvectors of both the conservative and dissipative operators) 


1. Spectrum of Ao- 

Let (To be the spectrum of Aq- We can split ao as follows: 

111 2 
(To — (Tq U (Tg, 


where 


cr, 


0 } i&Io 1 
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and Iq is a finite set, the multiplicity of is ruifl and is finite. 

^0 = {^k^}j=l,2,lkl>ko, 

and the multiplicity of = 1,2) is one. 

2. Eigenvectors of ^o- 

For eachi G Iq, we will denote by (j)\,l = 0, —1, a system of independent eigenvectors 

associated with G Ug. 

For each k G Z, |/c| > /cg, we will denote by = 1,2) an associated eigenvector of 
= 1,2) G erg. 

Moreover, since ^g is skew-adjoint, the system 

^0 = U {0i}|fc|>fcoj=l,2 

can be chosen such that jpQ forms an orthonormal basis of Ft. 


3. Spectrum of A. 

Similarly, let a be the spectrum of A. We can split a as follows: 

a = U (T^, 


where 

and I is a finite set, the algebraic multiplicity of Ki is mi and is finite, the geometric 
multiplicity is n*, with 1 < n* < m*. 

= {K}j=i,2,\k\>ko, 

and the multiplicity of = 1,2) is one. 


4 . Generalized eigenvectors of A. 

For each i E I, we will denote by k = 1..., n*, a system of independent generalized 

eigenvectors associated with Ki ^ W, which forms Jordan chains, i.e Sik > I, 
k = 1 , ...,ni, YlkU ^ik = ^i, 

{A - Kil)i)\k = ^ = 1, •••, ^ik, 

where we assume that = 0. 

For each k G Z, |A;| > /cg, we will denote by fiKj = 1,2) an associated eigenvector of 
XiU = 1,2) G 
The system 

^ ~ U {'4^k}\k\>ko,j=l,2 

is chosen such that any ^>1 G J^, \k\ > ko,j = 1,2, satisfies = 1- 
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3.2 Eigenvalues of A. 

Let A G C* and t/ 7 ^ 0, ?7 = (m, v , y, z, r], 7 ) G D{A) such that 


(20) AU = XU. 

Then 7 = n(l), 7 = z{l) and (n, v, y, z) is solution of 


v(x) = Au(x), X G (0; 1), 

Uxx(A + yx(x) = Av(x),x G (0; 1), 

z(x) = Ay(x),x G (0;1), 

yxx(A - bux(x) - by(x) = Az(x),x G (0; 1), 

m(0) = 0, 

2 /( 0 ) = 0 , 

A^m( 1 ) + ki(ux(l) + 2 /( 1 )) + k2Au(l) = 0 , 

, ^^ 2 /( 1 ) + hyxil) + k4Ay(l) = 0 . 


Eliminating v and 2 ; implies that solving (20) is equivalent to solving: 


' (Uxx + yx-A‘^u)(x) = 0,x e (0;1), 

(**) (Pxx - bUx -by- A^y){x) = 0,x G (0;1), 
. . {in) m(0 ) = 0, 

^ ’ 1 {iv) 2 /( 0 ) = 0 , 

(n) A‘^u{l) + ki{ux{l) + 2/(1)) + k 2 Au{l) = 0, 

, (m) A‘^y{l) + ksyxil) + k^Ayil) = 0. 

From {{) and {ii), it follows that u is solution of 


( 22 ) ~ 2A^M)f^(a;) + (A^ + h)A^u{x) = 0 

(cf. ([I^ with a = 1 and A replaced by (—iA)). 


Denoting by fi, t 2 , ^3 and the solutions of the characteristic equation 2A^r^+A^(A^+6) = 0, 

i.e. 

(23) ti{A) = ti = AA\/ iVb + A,t2 = 23 (A) = h = AA\J -lAb + A, 24 = - 23 , 

the general solution of ( 2 ) and {ii) is proved to be given by 

4 4 

(24) u{x) = ^ y{x) = ^ CidA'"', 

i=l i=l 


where Cj G C, i = 1, ...4 and 
(25) 


-v + «? v-i= -y + ti 

CLi — -;-, (I2 — -;-, U3 — -7-, 0,4 — 


2l 


2 l 


23 


h 


The values for dj, i = 1,..., 4 come from (2), using the expression for u given by (24). 
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Note that (20) and © imply 3fJ(A) < 0. In the proof of Theorem |2.5| the absence of purely 
imaginary eigenvalues is proved. Thus 3fJ(A) < 0 and ti does not vanish nor The coefficients 
di, ^ 2 ; dz ci-nd are well defined. 


Therefore the boundary conditions {Hi) — {vi) are equivalent to the system 


where 

(26) 

(27) 

(28) 


/ 1 1 

9i{ti) gi{t2) 

X^e^^g2{ti) X^e^^g2{t2) 
\ A^e*i^3(fi) X^e^^g^{t2) 


1 1 \ 


^ Cl ^ 

CO 


C2 

X^e^^g2{h) X‘^A*g2{tA) 


C3 

A^e*3^3(f3) A2e*4^3(f4) / 


V C4 / 


0 , 


9i{i) 

92{t) 

9z{t) 


A2 

■‘ + T' 

^2^ + {ki + t)\ 

Xt ’ 

{-e + A2) {U + A(fc4 + A)) 

XH 


Multiplying the third and fourth lines of the previous system by 


1 

A2 


, this one is equivalent to 


(29) 


( ^ ^ 

9i{ti) 9i{t2) 

e*i^2(ti) e*2^2(^2) 

\ e^^9z{ti) e^^g^{t2) 


1 1 \ 


/ Cl ^ 

gi{tz) giiti) 


C 2 

e^^g2{h) e^^g2{U) 


C 3 

e^^gzih) e^^gsiU) y/ 


V C4 / 


0 . 


Let M(A) be the matrix of the previous system and C = ( 01 , 02 , 03 , 04 )*, then we deduce that 
A G C (5R(A) < 0) is an eigenvalue of A if and only if A is solution of the characteristic equation 


(30) det(M(A)) = 0 -v^/(A) = 0 , with /(A) :=--det(M(A)). 

166 

(The division by (—166) simplifies the expressions calculated in next subsection for the asymp¬ 
totic analysis.) 

If A is an eigenvalue of M, an associated eigenvector has the form 


U = (u. Am, y, Xy, Am(1), A|/(1)), 


and is given by C a nontrivial solution of (29) and formulas (24)-(25). Moreover the geometric 
multiplicity of A is equal to the dimension of the kernel of M(A). 


Note that the expressions of g 2 and g^ depend on the values of /c 2 and k^. Thus the eigen¬ 
values and eigenvectors of Aq are different from those of A. 


3.3 Asymptotic analysis 

In this part we study the asymptotic behaviour of the large eigenvalues which are proved to he 
in the strip 

B = {XeC:-a< 3fJ(A) < 0}, 
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where a > 0 is fixed and chosen large enongh. 

The large eigenvalnes are also proved to be simple and the asymptotic expansions (36) and (37) 
are established. 

We first start by: 

Lemma 3.2. (Asymptotic behaviour of the characteristic equation) 

There exists a > 0 such that the eigenvalues of A are in the strip 

B = {XeC:-a< 3fJ(A) < 0}. 

Moreover the characteristic equation admits the following expansion 


(31) 




where /j, i = 0,3 is a bounded function on B given by (35) below 


Proof. Firsh if A is an eigenvalne of the operator A associated to the normalized eigenvector 
U, from (11), 0 > 3fJ(A) = —> —^2 — ^ 4 , since ^ ■ bP and ^ ■ hP are both 

_ ki % ki Kq 

smaller than Ilt7lb = 1- Hence the existence of a. 


Fnrthermore i = 1..., 4 is bonnded as |A| —> oo, where tj = ti{\),i = 1,4 is given by (23). 
By Taylor series it holds 


(32) 


iy/b b ib2 „ , * 

i. = A+ —+ —+01^^ 


(33) 


i\/b b ib^ ^ / 1 

is = A- —+ - + ^+ 0 ( 3 ^ 


Inserting (32) and (33) into (29) and nsing Taylor series, after long calculations we get 


M(A) = M(A) + 0(^-J, 

where M(A) is a matrix which only contains terms of order 1, — or —. Computing the determi- 

A A^ 1 

nant of M(A) and keeping only the terms of order less than or equal to —tt, we get after lengthy 
calculations 


(34) 




where fi,i = 0, ...,3 is a bounded function given by 
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/o(A) = 

/1(A) = — — (2(^2 + k^) — (/ci + ^2 + ^3 + k/i) + e *^(^1 + ^3 ~ ^2 ~ ^ 4 )) ) 

/ 2 (A) = — — {—4(6 + 2A;i/c 3 — 2fc2^4) 

lo 

+ (36 — 4fciA:3 — Ak2kz — 4,kiki — Ak2k/C)e^^~^^^ 

+ (36 — 4:kiks + Ak2k3 + Akik4^ — 4fc2^4)e“*^“*^ 

+ (—6 + 2iVbki — 2iVhk'i)e^^~^^ 

+ {—h — 2iVhki + 2iVhk'i)e~^'^^^^} ^ 

/ 3 (A) = -^{-ih(h + ki) 

lb 

+ ^6(7/ci + 6^2 + 3^3 + 

+ -^6(7fci - 6A;2 + 3 /C 3 - 6A;4)e“*^“‘^ 

+ {—bk2 — 2iVbk2k3 — 6^4 + 2iy/hkiki)e*'^~^^ 

+ (— 6 / c 2 + 2iVbk2k^ — bk^ — 2iVbkik^e~^'^^^^}. 


□ 


Lemma 3.3. (Asymptotic behaviour of the large eigenvalues of A) 

The large eigenvalues of A can be split into two families (A'^,)^^^ |fc|>fco' ^ ~ ^ chosen 

large enough.) The following asymptotic expansions hold: 


(35) X\ = ikn + o{T), \\ = ikn + o{T). 

Either = A^ and this root is of order 2, or 7 ^ A^ and these two roots are simple. 


Proof. The multiplicity of the roots of /q given by (35) is two and A is a root of /q if and only 
if 


3k G Z, (6i + t3)(A) = 2ikTr. 


4A 


Since (64 + 63 )(A) = 2A + — + o I — , we deduce that, for each A; G Z, with \k\ large enough, 


corresponds a double root of /o, denoted by A^ which satisfies 


Afc = iA;7r + 0(-). 

We will now use Rouche’s theorem. Let Bk = B{ik7r, rk) be the ball of centrum ikn and radius 
1 

Tk = 


kdP 


and A G dBy^ (i.e A = i/cvr + r^e*®, Q G [0, 27r]). Then we successively have: 
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(^1 + t3)W 


= 2ifc7r + 2rfce*® + O(-), 


k' 


g(4l+t3)(A) ^ ^2r,e“»+0(^) 

= 1 + 2rfce*^ + O(r^), 


and 


/o(A) = (1/4)(1 - 2rfce*^ + 0(r^))(2rfce*^ + 0(r^))2 
= (1/4)(1 - 2rfce*'^ + 0(r^))(4r^e^*® + 0(r|)) 

= + O(r^). 

It follows that there exists a positive constant c such that 

VA e dBk, |/o(A)| > crl = 


Then we deduce from (31) that |/(A) — /o(A)| = 0(y) = follows that, for \k\ large 

A K 

enough 

VAeE,, |/(A)-/o(A)|<|/o(A)|, 

hence we get the result. □ 

Remark 3.4. Since the imaginary axis is an asymptote for the spectrum of A, then System 


(29) is not uniformly stable. 


Remark 3.5. Obviously the previous asymptotic analysis of the spectrum is not necessary to 


deduce that System (29) is not uniformly stable. Indeed, using the compact perturbation result 


of Russell (see m), we directly see that the dissipative system n2W is not uniformly stable. 

More information concerning the asymptotic behavior of the spectrum of A is given by: 

Proposition 3.6. (Asymptotic expansions for the eigenvalues of A and Ao) 

Assume Condition 

(Cl) : ki k^ or Vb ^ 2k'K, k G N*. 

Then the large eigenvalues of the dissipative operator A are simple and can be split into two 


families (A;(,) 


k^X,\k\>ko 


,j = 1,2, (ko G N, chosen large enough.) 


Moreover, we have the following asymptotic expansions for the eigenvalues of A: 


Xl=tk7r + t^-^ + o(^), 


^.02 /^2 
'T 


At = + »— - ^ + “ I - 


(36) 

(37) 

where aj G M, (3j > 0,j = 1,2. 

If Condition (Ci) above is still assumed, the large eigenvalues of the conservative operator 
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^0 are simple and can he split into two families \k\>k° (^o ^ c/iosen large 

enough) with the following asymptotic expansions: 


(38) 

(39) 


with the same aj as above. 
(cf. Figure 1 of Section]^) 


Xl^ = ^klT + + O ) , 


Xlfi = ^klT + + o ) , 


Remark 3.7. The explicit values for ai and a 2 are given by (47), (45) and (46). They only 


Ulo 


depend on the values ofb, ki and k^. As for fdj, it is defined by fdj := j = 1; 2, with cj{ and 




UJ 2 given by (49) and (50) 


Proof. Step 1. 

Let Afc = Xl, with j = 1 or j = 2. From (36), it follows Xk = ikn + e*,, where Ck = o(l). 
Using (32) and (33) leads to: 


U + ts = 2zfc7r + 2efc 


ib 


+ o{e,) + o(jPj+o (I) , 


which implies: 


(40) 

(41) 


gtl+t3 ^ X 


ib 6^ ibck 


Akn 32k‘^7i‘^ 2k7r 

= 1 , A _ 

4A;7r 32A;^7r^ 2k7r 


+ 2e, + ofe) + o (^) + o (j) , 
- 2t, + o(t») + o (^1 j + o (I) , 


Similarly it holds 


and we deduce that 


. rr ib^!"^ ( 1 


(42) 

(43) 


= e*''® + 


g-ti+ts _ 


z63/2g*\/l> 

8A)%2 

^^3/2g4v4) 

8A:%2 +""1^ 


Using (31), inserting (40)-(43) into fifXk) and keeping only the terms greater than or equal to 
0( —), we obtain after calculations 
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( 44 ) /(Afc) = el-i-fij- 72 ^ + o + o{el) + o{j) = 0 , 

where 


( 45 ) 


7 i = 


b + 4 (/ci + ks) 
Att 


( 46 ) 72 


— 8 b + 6^ + 8 bki + 86/^3 + GAkiks + 8 bcos{\/b) + lQ\/b{ki — ks) sin(-\/6) 

647r2 


Multiplying ( 44 ) by k^ leads to: 


{kekf - i'liikek) - 72 + o(l) + o{kek) + o{Pel) = 0. 
Thus kek is bounded and 


{kekY - hi{kek) - 72 + o(l) = 0 . 

The previous equation has two solutions 

= ^(71 - \/ 7 ? -472) + 0(1) or kek = ^(71 + \l ll - ^12) + o(l)- 

Denoting by 


(47) 

it holds: 


7i - 7/7? - 472 , 7i + 7/7? - 472 

«! =- - - and 0(2 = --, 


. cq 


k 


. 0:2 


efc = + o r or ek = %-r + o - 


k) " fc 

Note that, if Condition (Ci) holds, ai and 0:2 are real numbers and ai 7^ a2- Indeed 71 G 
and it holds 

2 6 + 2(/ci — ^3)^ — 6 cos(7/&) — 2a/6(A:i — fcs) sin(A/6) 

7i - 472 = 


27r2 


1 


/ -| \ ‘^ h 

2 (ki — ks — -y/bsm{y/b) j + ^ (^{cos{\/b)Y — 4(cos(a/ 6)) + 3) j 


> 0 


for all ki > 0 , k^ > 0 ,b > 0 . Thus 7^ — 472 = 0 , if and only if ki = k^, and a /6 = 2 kn, fc G N*. 


Now it must be proved that near ifcvr, there are exactly two distinct roots, for \k\ great enough. 

For that purpose we consider Tfc tl 
and the polynomial pk defined by 


1 — C^2 

For that purpose we consider F^; the disk of center = ikir + i—- and radius Vk =-^- 

k 2 k 
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Pk{z) = {z- ikirf - i 7 i 


2 ; — ikn 
k 


72 


k^' 


The roots ofpfc are 2 ;° and ikn+i— (it holds af—jiai+'y 2 = 0 since ai+a 2 = 71 and aia 2 = 72 )- 

K 

Gl 1 

But ikir+i— does not belong to T^, if \k\ is large enough. Let ;2 = zl + — {02 — ai)e^^, 9 G [0, 2n] 

K ZdK 

any element of dVk- Then pk{z) is proved to be: 

(e*^ - 2) (7^ - 472) 


Pk{z) = 


4k‘^ 


thus there exists a positive constant c independent of k such that 


\Pk{z)\ > p, V2; e dVk. 


On the other hand, using (44) we get 1 /( 2 ;) — Pfc( 2 :)| = o ( ^ ) • Therefore, Rouche’s theorem 


implies that / has only one root in Tfc, if k is large enough. Finally, we have proved that 
the large eigenvalues of A are simple and can be split into two families with the following 
expansions: 

Al = zk7r + t^ + o , Al = zk7r + z^ + o . 

Note that the eigenvalues of the conservative operator Aq have the same asymptotic expansions, 
since ai and 0:2 are independent of the values of k 2 and ^ 4 . 


Step 2. 

Since for j = 1, 2, aj G M, we need one more term in the expansion of A^A = 1, 2. 
From Step 1, the expansion for j = 1 or j = 2 is: 


= + + P 


where ei = o(l) 
Using I 
after a 


Using (31), Taylor series and simplification in the term of order — coming from Step 1, we get 

k^ 


Tong calculation 


(48) 




/(z/c7r + ^^ + ^) = ^{uiei + {eif)+u;i^ + o(^ ) =0 


where uj G zM, Z = 1, 2 and is given by 


(49) 


/ A — 

(Ui — 


i{b + Aki + 4^3 — Soj-tt) 

47r 


= Th/ 7 i -472, 


(50) 


^2 = P^(73 - 87r(A;ifc2 + k^k^aj), 
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j = 1 , 2 , where 


(51) 


73 = h{kik2 + ^ 3 ^ 4 ) + 8 / 01 / 23(^2 + k 4 ) + 2Vb{kik2 — / 03 / 04 ) sin {Vb). 


Since we assume (Ci) then uj ^ 0 (see the remark just below (46)) and we deduce from (37) 




J I 


uf k 


uj-. 


that e], =- 17 + *^( 7 ) - Setting I3j = then it holds jSj G M and (36) holds. Since all the 


Wi 


eigenvalues of A are on the left of the imaginary axis, necessarily (5j > 


IT 


Note that, if /02 = /04 = 0 (associated conservative operator Tlo), 73 = 0 and thus, 0 J 2 and 
f3j vanish as well. 


Now, if (/c 2 , / 04 ) 7 ^ (0, 0) (dissipative operator A), (3j 7 ^ 0, j = 1, 2 as it is proved below. 


Step 3. 

Assume that (/ 02 , / 04 ) 7 ^ (0, 0) and (3j = 0,j = 1, 2. Then = 0, j = 1, 2, thus 

a ■ = _ — _ 7 = 12 . 

^ 8TT{kik2 + ksk^)' 

But, since ai + 0:2 = 71 and ai • 02 = 72, it holds: 

-7i«i +72 = 0, j = 1,2. 

It follows 


73 


-7i; 


73 


+ 72 = 0. 


87r(/ci/o2 + /C 3 / 04 ) / 87r(/ci/c2 + / 03 /C 4 ) 

We multiply the previous identity by 16(/ci/o2 + k^.k^Yn'^ and use (46) and (51) to get: 


74 + 75&cos(\/ 6) + 76&sin^(v^) + 777^ sin {\fh) = 0 , 

where 74 = — 2 b{kik 2 + k^kY)"^ — 16/oi/c2/o3/o4(/oi — /c 3)^,75 = 2 (/ci/o 2 + k^k^Y, 
76 = (/01/02 - ksk^Y and 77 = 16/oi/o2(/oi - k^k^k^^. 


Now, using the fact that 74 — 472 > 0 or equivalently 27r^(7^ — 472 ) > 0 (which is true if 
and only if Condition (Ci)), it holds 

b + 2 kl — 4/04/03 + 2/03 — 6cos(-s/6) — 2 Vb{ki — /03) sm {'/ b ) > 0. 

Thus, using the dehnition of 77, 

77*7^ sin (a/6 ) < 8/04/02/03/04(6 + 2(/oi — k^Y — 6cos(\/6)). 

We get after simplihcations 
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0 


= 74 + 75^cos(v^) + 76fesin^(A/6) + 

< {2(kik2 + — 8kik2k3k4){—b + bcos{\/b) + {kik2 — /c3fc4)^6sin^(-\/6)) 

< b{kik 2 — kzk^‘^{—2 + 2 cos(\/6) + sin^(A/6)) 

= -Ab{kik2 - k^ikiY ■ 

Since this inequality never holds, the assumption f3j = 0, j = 1; 2 does not hold either. 

□ 

Now, if Condition (Ci) does not hold, the calculations are different (and long). The details are 
not given here. The results are given without proofs. 

Proposition 3 . 8 . (Asymptotic expansions for the eigenvalues of A and Aq - particular cases) 


1. Case ki = /ca, k 2 Y ^4 ^ p G N*. 

The large eigenvalues of the dissipative operator A are simple and can he split into two 
families {^i)ke:Z\k\>ko A ~ (ko G N, chosen large enough.) Moreover they satisfy the 
following asymptotic expansions: 


(52) 

(53) 




ikiT + 


i(2ki + p^TT^) 
2k'K 


ikn + 


i{2ki + p^TT^) 
2k'K 


kik2 / 1 \ 

hh / 1 A 


(cf. the table and Figure 2 of Section 


2 . Case ki = ks, ^2 = ^4 7^ 0, b = dp^vr^, p G N*. 

The large eigenvalues of the dissipative operator A can he split into two families {^k)ke:Z \k\>k ’ 
j = 1, 2, (ko G N, chosen large enough.) Moreover they satisfy the following asymptotic 
expansions: 


A{ = ikiT + 


i {2ki + p^TT^) 
2kTT 


kik2 ijasj - 2Akikl) / 

^ 2^2 24:k^7r^ \k^) ' 


where G M,j = 1,2 are given below. 


3. Case ki = k^, k 2 = k 4 = 0, b = dp^vr^, p G M*. 

The large eigenvalues of the conservative operator Ao can be split into two families 
k&\k\>k° A = I 52 , (k^ G N, chosen large enough) with the following asymptotic 
expansions: 


where G M,j 


A-^’° = ikii + 




1 , 2 are given below 


2k'K 2d/c^7r^ \k^) ' 
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03,1 = —2Ak‘l — Ski ~ SGkip^TT^ + 9p^7r'^ — 12p7r-JAkf — 43A;fp^7r^ — Akip'^ir'^ + p^n^, 


03,2 = —2Akl — Ski — 36fcip^7r^ + + 12p'n\J Akf — ASklp'^n'^ — Akip^n"^ + p®7r®. 

Note that, if Akf — 43/cfp^7r^ — Akip'^n'^ + p^tt^ ^ 0 then A^’° ^ A^’° for k large enough. 
Idem for Af and A^ of the previous case. 


4 Riesz basis 


In this section, it is proved that the system T of generalized eigenvectors of the dissipative 
operator A (introduced in Theorem 3.1) forms a Riesz basis of H. To this end, we use Theorem 
1.2.10 of [I] which is a rewriting of Guo’s version of Bari Theorem with another proof (see I3)- 
For the sake of completeness. Theorem 1.2.10 of jT] is recalled : 


Theorem 4.1. Let A he a densely defined operator in a Hilbert space H with compact resolvent. 
Let {(fn}'i^=i he a Riesz basis of H. If there are two integers Ni, N 2 > 0 and a sequence of 
generalized eigenvectors {'fn}ff=Nx+i of A such that 

00 

^ ^ 11 0n+Af2 'fn+Ni 1 1 2 ^ ) 

n=l 

then the set of generalized eigenvectors (or root vectors) of A, {'ifn}f(Li forms a Riesz basis of 

n. 


The family Hq of eigenvectors of the conservative operator Alo is an orthornormal basis of the 
Hilbert space H. Thus, it is enough to show that the eigenfunctions of Ao associated to the 
eigenvalues Af^ E (Tq and those of the dissipative operator A associated to the eigenvalues 
Af. E cr^ are quadratically close to one another. 

Theorem 4.2. (Riesz basis for the operator A) 

For any j E |1;2|, it holds: 

j: \\4-<r,ru<oo. 

\k\>ko 

Thus, the set IF of generalized eigenvectors of A forms a Riesz basis ofFi. 

Proof. Step 1. 

Since fji lies in Ft, it has six components (see Section]^. Let us write fjl := (n;^, vf, yf, zf, pf, 7 ^) 
and let us first prove that 

171 = 0 ( 1 ), 171 = 0 ( 1 ). 


(54) 

From 0 , it follows 
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(55) 


3? < >H= -^Wk? - ^l7fcl 


Now, 3? < >H is also equal to 3fJ(A;i) = + o 


(37). Hence (54). 


O I — 1 due to 


(36) and 


Step 2. Projection. Let j = 1,2 and k, \k\ > ko be fixed and denote by the or¬ 
thogonal projection on the orthogonal space of the 1-dimensional space directed by (jy’p.. 

Clearly there exists «;(, which can be supposed to satisfy > 0 without loss of generality, such 
that 

(56) V'fc = + Rl 


where Rl = PH'il’i) and Rl is orthogonal to 0;(,. 
Thus, due to Lemma [4.3| given later. 


Then, 3c(,, real number bounded with respect to k, such that 




Step 3: {(l>i}\k\>ko {'4’k}\k\>ko quadratically close to one another. 

Using Step 2, 


Hence 






|A :|>0 




111 2 
kWn 


< oo. 


□ 


4.2 


Lemma 4.3. (Technical Lemma for the proof of Theorem 4-2 ) 
The vector Rj,, defined in Step 2 of the proof of Theorem 
J = l,2. 

Proof. Using (|56|, it holds, for j = 1; 2: 


satisfies: \\Ri\\n = O(-), for 


(.Ao - — (Tlo - + Ri) - ai{Ao - + (A - K)iRi) 

= k 2 r]le 5 + ^47^66 (this follows from (9)). 

Since Aq and Pi commute, then applying Pi to the previous identity, we get 


(A - Xi){Rl) = hviPlie,) + k.^iPlie,). 


19 












Thus 


Ri = hviiA-K) ^Plie5) + h'yi{Ao-Xi) ^P^(e6). 

Writing 65 in the orthonormal basis J^q, it follows 




P’M) = E 

+ [< eg, 0; >H ■<Pi+ < 65 , 0 ^^ >n >] 

\l\>ko,¥k 

+ < 65,0-^’'’^ ■01"''^) 


where the exponent j is dehned modulo 2. 

Then 

||(A-Ai)-‘(P^(e5))||„ < ^ 


< 65, (Pi>n\- 


iPi - K\ 


E 

|Z|>fco,i^fe 


< 65, 0f >-^ 


+ I < eg, 0; ^^ >-H 
+ I < eg, 0-^^^ 




> 


1 


E C'||eg||'H + I < eg, >'n 


and similarly 

WiAo - X{)~^{Pi{eG))\\n < C\\eG\\n + \ < eG,(j){'^^ >n 


lA^’^-Ail’ 


lA^’^-Ail 


The existence of the constant C independent of k in the latest expressions comes from the fact 
that I is different from k in the sum. Indeed the behaviour of |Af*^ — A'^l and that of |A/’*’^’° — A'^l 
are given by (36), (37), (38) and (39). They are both bounded from below by a constant inde¬ 
pendent of k and 1. Note that this still holds in the particular cases described by Proposition 

The expression ~ A{| is not bounded from below by a constant independent of k. The 

same asymptotic expansions prove that, for j = 1; 2: 


(57) 


|Ar'“-Ai:|=0(-). 


Thus, using (54), the result follows as soon as it has been proved, for j = 1; 2: 


(58) 


< eg,0i^^ 


>H I = 0{\) and I < e6,0fc^^ >'h \ = 0{\ 


'k' 


k' 


Since 0^ := (4’°^^°) 2/fc^ 4’°) 7fc°) ^ D{A), and eg = (0,0, 0, 0,1,0) then 


(59) 


< eg, 04^ >n= vf = 


and (4’°, 4’°) 4’°) 4’°) solution of System (21) with A = A^.’^ = *4’°, 4’° G M. In particular 
(i) is 
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+ {yf)^ = -{hf)'- 


u 


k ■ 


For simplicity, the indices and exponents are dropped from now on. 


{u^x + yx){x) ■ {u^ + y){x)dx =-/ h‘^u{x) ■ {Ux + y){x)dx 


= -h^ Q u{x)ux{x)dx + J u{x)y{x)dx'^ . 


Integrating by parts, it follows 


u{x)ux{x)dx = — Ux{x)u{x)dx + |m(1)|^ + \u 


and, due to {iii) of System (21): 


23 ? < / u{x) ■ Ux{x)dx\ = \u{l)\ 


And thus 


(60) {uxx + yx){x){ux + y){x)dx^ =-h^ + 2^\ I u{x)y{x)dx 

On the other hand, after an integration by parts, it holds: 


H“ Vx) ' "b y^dx J {^Ux “b 2 /) ‘ {^Uxx “b yx^dx 

+ |mx( 1)+2/(l)P - |mx( 0)+2/(0)p 


which implies 


(61) {ux:x + yx){x)-{ux + y){x)dxj = \u^{l)+y{lX-\ux{0) + y{0)\^. 

Now, using V = Xu and z = Xy (cf. the system just before System ([^), ([^ and (61) imply 


(62) 


Ma:(l) + |/(1)| + h |m(1)| = |mj;(0) + ?/(0)| - 23? < / v{x)z{x)dx 


Then, (n) of System (21) with k 2 = 0 {(pl is an eigenfunction of ^o) leads to 

kj ■ |M:r(l) + 2/(l)p = 


And 

(63) 


A|‘‘.|u(i)p = ;i-‘.|u(i)p = ft2.|t,(i)p 

«,(!) +s(l)P + ft"Wl)P = (l + |)|!.(l)p. 


Using (iv) of System (21) as well as the trace Theorem applied to Ux implies that there exists 
a constant Ci such that: 


(64) 


|mx(0) + 2/(0)p = |mx(0)P < Cl (^WuxxWh^n) + \\u. 
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Now (i) of system (21) gives: 


Ur 


~ 11'^^'“ y^Wl^in) — l-^l^ ■ ll'*^lli2(o) + lba^lli2(o) 

< + H(l>\\n < UWn + b ■ U\\n <l + b. 


And 


— ll'^a; + 2/|li2(n) + ll2/lli2(Q) 

1 b 

< HWh + j^\\^\\h{n) < 1 + |yj2 ^ 


- /12- 


Using successively the two previous estimates in (64), the Cauchy-Schwarz inequality applied 


to the last term of the right-hand side of (|62|), (|63|) and (|59|), we get the hrst result of (58): 

1/2 


(65) 


< 65, 0i >n I < 


1 ^ 1(2 + p-) + 2 


fe2. 

fc2 


Indeed, by dehnition, /i{’° is the imaginary part of A{’^ which behaves like k for large values of 


k (cf. Propositions 3.6 and 3.8). 


To end this proof, let us give the sketch of the proof of the second estimate of (58). The 


ideas are similar to those developed just before. That is why the details are not given here. 
It holds < ee, >h= 7fc° = with the same notation as before. 

Integrations by parts allow to write the analogous of (62): 


(66) |i/a;(l) - 5 m(1)| + {h -6)|i/(l)| = |i/3,(0) - &m(0)| + 2b{h - 5) ■ 3? <j j yudx 


Long calculations, using System (21), the Cauchy-Schwarz inequality as well as the trace The¬ 


orem applied to i/x, lead to the existence of a constant C 2 such that: 

^ + 1 - |^(1)P < b{{b + 2)C'2 + 26) 6^(C2(6 + 1) + ‘^b)— + —|u(l)p. 


Using (65), it follows: | < e6,0|, >-^ | = O(-). 

fC 


□ 


5 Polynomial decay rate of the energy 


The energy is already known to be not uniformly stable (cf. Lemmas 3.2 and 3.3 and the 
remarks just below the lemmas). It is now proved to decay polynomially. To this end, the 
solution is explicitly expressed using the Riesz basis T of generalized eigenvectors of A (cf. 


Theorem 4.2). 


Theorem 5.1. (Polynomial demy rate of the energy) 

Assume that a = 1 in System Then there exists a constant C > 0 such that for any 

initial datum Uq G D{A), the energy of the system rewritten as (10) satisfies the following 
estimate: 
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E{t) < C ■ 


t 


where E{t) = - \\U{t)\\l^. 


Proof. Using the Riesz basis T (cf. Theorem 4.2), we can write 


^ik 


^0 = Y1 (^oyi^r 

\l\>ko,j=l;2 


The solntion of (10) is: 


u(t) = Y.^ 

iei 


Kit 


iGl k=l 1=1 


E E (E(“») 

/c=l 1=1 \p=l 


tp 


-I 


{p-l)\ 




+ Y1 

\l\>ko,j=l;2 


Since is a Riesz basis, there exists a positive constant K snch that the energy satishes, for 
any f > 0; 


E{t) < K 

-maxlP”*-!;!} ■ 

( A 

EEiwtP 

+ 


i&I 

\fc=i 1=1 / 

\l\>ko,j=l-,2 


Using the asymptotic analysis performed in Propositions 
all |/| > fco,j = 1,2, it follows that 


3.6 


and 


3.^ 


and since 3fJ(A;^) < 0, for 




E{t) < R'maxig/lf™* 1} • e 2mmig/(»(|Ki "‘EEEi(“»)i'ti 


E 




i&I k=l 1=1 


\l\>ko,j=l-,2 


P 




where /3j,j = 1,2 are positive constants. 

Now, if /d > 0 is hxed, the fnnction u ■ is a bonnded function on M’*'. And 

Til ^ik 

EEEllf^olli. “'d E P■ 

iel k=l 1=1 \l\>ko,j=l;2 

Hence the result. 


□ 


6 Numerical validation 


The asymptotic behavior of A^, given by Propositions |3.6| and |3.8[ can be numerically validated. 
For instance in the case ki = = 2, k 2 = 1, k^ = 5, b = dvr^, (hrst case of Proposition 3.8), 

we have calculated numerically some large eigenvalues near the imaginary axis. From (52) and 
( [5^ it holds, in that case: (3j, j = 1,2, with 
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-1.01321. 


f3i = -^ 


TT^ 


-0.202642, 


/S2 = -^ 


TT^ 


The table below confirms this behavior. 


k 

200 

400 

600 

800 

1000 

emi 

-0.202667 

-0.20265 

-0.202652 

-0.202588 

-0.202637 

k'^^RXl 

-1.01303 

-1.01317 

-1.0132 

-1.01324 

-1.01314 


The figures hereafter represent the eigenvalues of A in two cases: 


Proposition 3.6 and Figure 2 to the first case of Proposition 3.8 


Figure 1 corresponds to 


•• 
• •' 

•Jl- 

• • 




• • 


Figure 1: Eigenvalues of ^ if a = 1, 6 = 2, fci = 1, ^2 = 2, A:3 = 3, ^4 = 2 



• « 
• • 


• • 
• « 
• • 


Figure 2: Eigenvalues of ^ if a = 1, 6 = ki = = 2, ^2 = 1, k^ = 5. 
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